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Extremely strong magnetic fields change the vacuum index of refraction. This 
induces a lensing effect that is not unlike the lensing phenomenon in strong gravita- 
tional fields. The main difference between the two is the polarization dependency of 
the magnetic lensing, a behaviour that induces a handful of interesting effects. The 
main prediction is that the thermal emission of neutron stars with extremely strong 
magnetic fields is polarized - up to a few percent for the largest fields known. This 
potentially allows a direct method for measuring their magnetic fields. 
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1 INTRODUCTION 

The most intense magnetic fields in the universe are found 
near neutron stars. S o on after the d iscovery of radio pulsars 
(Hewish et al. 1968), Gold (1968) argued that the steady 



increas( i in their periods could be explained it they were 
strongly magnetized neutron stars with dipolar fields oi 
~ 10 12 G. Since that time, the variety of radio pulsars has 
expanded. These neutron stars have magnetic fields rang- 



panions have been found near these objects. Although no 
single picture has yet emerged to explain AXPs, a strong 
candidate model is that a neu tron star cooling throug h a 
strongly magnetized envelope ( Heyl fc Hernquist 1997q) or 
the decay of a neutron star's strong magnet ic field ( |Heyl & 



Kulkarni 199£ ; Thompson fc Duncan 1996 ) powers the X 



ray emission from these objects. Furthermore, such a strong 
magnetic field (B ~ 10 15 G) can also acco unt for the ob- 



ing frorh 10 G lor th e millisecond pulsars to 2 .1 X ID G 



served spin -down rate of these objects (e.g. Vasisht & Got 



for PSR 0157+6212 (Arzoumanian et al. 1994), the most 



thelf 1997) 



strongly magnetized radio pulsar discovered to date. 

A dozen years after the discovery of radio pulsars, the 
increased scrutiny of the sky at 7-ray energies uncovered a 
new class of objects, the soft-gamma repeaters (SGRs). In 



Intense magnetic fields of B 10 12 G strongly affect 
physical process in and near neu t ron stars - ranging f rom 



cooling ( Heyl fc Hernquist 1997q; Shibanov et al. 1995), to 



t EuccoEcful modol for thoso onorgotic eourcos, 



the mo 

Crust r>F a sfcEQIlgljZ magnptw.pd (R 



H) npntrmi 



star or "magnetar" periodically fractures, releasing a large 
amount of energy into th e magnetosphere which is ob served 
as soft-gamma radiation (Thompson & Duncan 1995). Until 



atmospheric emiss ion (Etajagopal, Romani fc Miller 1997 



^ Pavlov et al. 1994) and to the insulation of th e core (Heyl & 



Hernquist 1998b; Heyl & Hernquist 1998a; Schaaf 1990a 



Hernquist 1985). Yet stronger fields such as those asso- 
AXPs and SGRs alter the propagation of 
magnetosphere by way of quantum- 



the past year, the evidence for the intense magnetic fields 



ding SGRs was circumstantial , As- Gold (1968) 



gued fob radio pulsars^ Knnvplintnn pt. al (1 QQSa)| prprispjjr 



dated with 
light through the 
electrodynamic (QED) processes and may furthe r process 
the emergent radi ation (Heyl fc Hernquist 1997a|; Heyl fc 



measured the spin down of SGR 1806-20 and infer red a 
dipole field of 8 x 10 14 G. Kouveliotou et al. 



Adler 1971Q 



Hernquist 1997b; Baring 1995; Baring & Harding 199 



(1998b) 



sured P of SGR 1900+14 and find it to be consistent with 
B = 2 - 8 x 10 14 G. 

In the past decade, more sensitive X-ray telescopes have 
also discovered a new class of sources which may have ex- 
tremely strong magnetic fields. The anomalous X-ray pulsars 
(AXP) are soft X-ray sources and they are often associated 
with supernova remnants. Unlike the majority of X-ray pul- 
sars, these objects have periods of several seconds which 
steadily and quickly increases. Furthermore, no binary com- 



Not only will QED alter the spectra from these objects, 
it also creates a lens surrounding the neutron star which 
magnifies and reimages the surface. In a sufficiently strong 
field, the index of refraction for photons with their magnetic 
fields directed perpendicularly to t he local magnetic field ca n 
be significantly larger than unity (Heyl & Hernquist 1997b). 
Furthermore, it is a strong function of both the strength of 
the magnetic field and its direction. In this paper, we will 
examine the effects of magnetic lensing on the observations 
of strongly magnetized neutron stars. 
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2 INDEX OF REFRACTION 



3 _L-RAY PROPAGATION 



In the presence of a strong external field, the vacuum re- 
acts, b ecoming; magnetized and polarized. The index of re- 



fractior , magnetic permeability, and dielectric constant of 
the vacuum are straightforward to c alculate using quantum 



electrodynamic one-loo p corrections ([Klein fc Nigam 1964a 
Klein fc Nigam 1964bj, ^rber 1966|, lAdler 1971L |Berestet 



skii, Lifshitz fc Pitaevskii 1982, Mielnieczuk, Lamm & Val- 



luri 1988). Heyl fc Hernquist (1997b) calculate the index 



of refraction as a function of field strength to the one loop 
order. The magnetic field will strongly bend light as it prop- 
agates only if n — 1 ~ 1. 

For B > Bqed ~ 4.4 x 10 13 G, the index of refraction 
becomes significantly larger than unity if the wave has its 
magnetic field polarized perpendicular to the local magnetic 
field (_L-mode)Q To one loop, the waves with their magnetic 
field directed in the plane containing the local magnetic field, 
B, and the wavenumber, k (||-mode), have nil — 1 < a/(6n) 
regardless of the strength of the field (a, the fine-structure 
constant is approximately 1/137.04). This second mode is 
essentially unaffected by lensing. 

In fields where n± — 1 ~ 1, Heyl & Hernquist (1997b) 
find that n± — 1 oc B with sufficient precision. Specifically, 
for B > Bqed/2 
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where £, — B /Bqed and 9 is the angle between B and k. 7 is 
Euler's constant and In A is related to the derivative of the 
zeta function: In ,4 = 1/12 — (-1) = 0.2487... . The linear 
approximation is accurate to 20 % for the value of n — 1 for 
fields stronger than 3.9 x 10 14 G where n- 1 = 2.9 x 10" 3 . For 
this same field strength, n — 1 for the parallel polarization 
is ten times smaller. For stronger fields this ratio increases. 
Thus, we set riu = 1. 

The function f(B) = (n± — 1)/ sin 2 9 contains the por- 
tion of the index of refraction which depends on the strength 
of the external field. The linear approximation is given by 
f(B) — a/(6n) x B/Bqed = B/Bq. The second loop correc- 
tions to the index of refraction are expected to be s maller by 
a factor of a re gardless of the strength of the field ( Dittrich 
fc Reuter lggj). 



Note that the convention used here __ 
Adlor (1971). |Bcrcstctskii. Lifshitz fc Pitaevskii (1982) 



consis tent with 
(§130) and 



Heyl & Hernquist (1997b) in which the _L and || modes describe 
the polarization directio n of the ray's magnetic field relative to 
the local magnetic field. Baring (1995)| and other authors use the 



In the limit of geometric optics, light ray s obey Hamil- 
ton's equa tions with H = u = cfc/n(k, x) (Landau & Lif- 



shitz 1987, §53). Since the index of refraction depends both 



on position and wavevector, the resulting equations are non- 
trivial, 



x = Vk.uJ and k = — V x lo. 



(2) 



Each of these gradients may be calculated for the _L-mode 
using the index of refraction, n± — 1 + [1— (k • B/kB) 2 ]f(B): 
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(4) 



As outlined in the previous section, the propagation depends 
on the polarization of the photon; thus, it is important to 
understand how the photon's polarization evolves as the di- 
rection of the field changes. 



3.1 Polarization evolution 

As the waves travel through the magnetosphere, they stay 
in the same polarization mode (parallel or perpendicular) as 
long as 



1+ni 



n± 



Ak 



< 1 



(5) 



where r n is the distance over which the magnetic field ro- 
tates by one radian and Ak is the difference between the 
wavenumbers of the two polarizations. Evidently, the adia- 
batic approx imation holds for the vicinity of the systems of 
interest (see Heyl fc Shaviv 1998 for further details). 

A photon's polarization will decouple from the magnet- 
ically induced polarization modes at a distance from the star 
which depends on the strength of the field at the stellar sur- 
face and the energy of the photon. Higher frequency photons 
decouple later; therefore, the observed polarization direction 
will depend on photon energy with higher frequency photons 
tracing the magnetic moment of the star projected onto the 
sky from later in the star's rotation. 

Generally, the wavelength of the photons is much 
shorter than the length over which the field geometry 
changes, so the decoupling occurs where the magnetic field 
is much less than the QED critical value. In this limit, the 
wave decouples from the field after traveling for a time, 



Couple « 4.6 ms^^j [w^G J lliTknJ ' (6) 

where B p is the strength of the magnetic field at the surface, 
R is the radius of the neutron star and E is the energy of 
the photon. 

The typical value of tdecoupio for high B neutron stars 
is such that decoupling will take place well within the light 
cylinder yet far enough from the stellar surface such that the 
local magnetic field where the photons decouple is aligned 
with the apparent magnetic axis, irrespective from where the 



6/5 



wave's electric field to denote the polarization. 
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photons originated. This simplifies the calculation of the two 
polarization images - The two principle polarization direc- 
tions at the observer are the || and _L modes even though 
the polarization directions from where the photons orgini- 
ated were different. 



4 THE APPEARANCE OF A MAGNETAR 

The nontrivial index of refraction depends on the direction 
and magnitude of the local magnetic field. Light rays will 
therefore travel in bent trajectories and the optical appear- 
ance of the magnetized neutron star will therefore not be 
trivial. From equivalent optical configurations (e.g. a black 
body embedded within a nonspherical medium of index n), 
it is clear that the images of the star will be distorted and 
can have a varying apparent surface area, even though it is 
spherical. 

We use the ray tracing algorithm to find the image of 
a magnetar having a dipole magnetic field. For each observ- 
ing direction characterized by an inclination angle above the 
magnetic equator, a different image may result. These im- 
ages can be compared to the unmagnetized case which ad- 
equately describes the || polarization. Each image can then 
be averaged in order to compose a light curve. And finally, 
all the images can be properly averaged together, to get the 
total emission properties of the NS. 

4.1 The ray tracing algorithm 

Using the Hamiltonian formalism and the related adia- 
batic and geometrical optics approximations (which hold 
extremely well in our limit), we know that each polariza- 
tion mode (either the |] or _L modes) will evolve separately 
according to eq. |^: 

= Vfec^j_ : || and k ± ,|| = -V x ^±,\\- (7) 

Since we are interested in the regime in which the index 
of refraction n« of the ||-mode can be assumed to be 1, its 
evolution becomes trivial. The image observed in the _L po- 
larization will not however be straight forward and its light 
ray trajectories will be bent. 

To construct the image observed by a distant observer, 
we can use the time reversal symmetry of the light trajectory 
and trace the light rays from the observer to the NS and not 
vice versa. We place a screen at a large distance from the star 
and divide it into pixels at the desired resolution (e.g. 300 2 ). 
From each pixel, we trace back the path of a light ray leaving 
in a direction perpendicular to the screen. The ray is then 
followed until it intersects the NS surface or until it is found 
to have an increasing radial component, implying that it has 
missed the star. The configuration is depicted in figure [l]. 

It is clear from Kirchhoff's law that if we look at a 
light ray trajectory that intersects the surface, then if the 
observer and the object are in thermal equilibrium, we will 
find that the specific intensity I (flux per unit sterad) of the 
ray going from the observer to the object should be equal 
to that of the opposite ray going from the star. Since the 
observer is at a large distance and has n — > 1 in its vicin- 
ity, the specific intensity is just I = B(T) — aT 4 /ir, where 
B(T) is the black-body source function (once the ray is in 
an n = 1 region) . Hence, if the emission from the NS is that 




Figure 1. The ray tracing algorithm. Ray trajectories are fol- 
lowed from a distant screen towards the star using the Hamilto- 
nian equations of motion. The rays can then either intersect the 
star or miss it. Since the magnetic field is assumed to be that of 
a dipole (aligned with the 2-axis), it has axial symmetry and dif- 
ferent images will be distinguishable by the observer's magnetic 
inclination angle i. 

of a black-body, then the rays coming from each element 
will simply carry a specific intensity of B(T) where T is the 
local effective photospheric temperature at the point where 
the ray intersects the surface. If the NS has a uniform tem- 
perature then the flux observed by an observer at a given 
direction is simply proportional to the solid angle subtended 
by the object. As we shall see, this apparent solid angle will 
be larger from all 47r possible observer directions implying 
that the total luminosity emitted at a given temperature 
is larger in the presence of a magnetic field. This does not 
necessarily imply that a magnetized NS will shine brighter. 
As we shall see in §5, the changed surface properties will re- 
duce the temperature to compensate for the higher emission 
efficiency, leaving an unchanged luminosity. 

4.2 Image properties 

A few resulting images can be seen in figure [| which depicts 
the magnetic latitude and longitude observed at several in- 
clinations for both a weak field and a strong one. We find 
that in the limit of a weak field (one for which the correction 
to the index of refraction is less than half), the effects of the 
magnetic lensing are to slightly distort the image through 
the imaging of otherwise inaccessible surfaces (surfaces that 
do not have a direct line of sight) and to have the images 
appear to subtend a larger surface area. In strong fields with 
corrections larger than a half to the index of refraction, on 
the other hand, the images are distorted enough as to have 
a topologically nontrivial image form in which some surfaces 
are seen twice (or more) with changing parities. 

The critical equatorial field needed to trigger the non- 
trivial behaviour is 2.8 x 10 16 G for which fo = /(-Bcquator) = 
0.25 on the equator and f(B po \ c ) — 0.5 on the poles. At the 
latter field strength (for which n± can reach 1.5), the group 
velocity is no longer a monotonic function of k, allowing 
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Figure 2. The magnetic latitude and longitude of the surface of the NS as seen with the _L polarization by a distant observer at different 
inclinations above the magnetic equator (with spacing of 10° in each direction). The equatorial magnetic field strengths of 2.3 X 10 16 G 
and 2.3 X 10 17 G corresponding to /o = 0.2 and /o = 2.0. An observer of a "weak" field NS will see a larger image and more than 2-7T 
sterad of the surface. If the NS is a black-body, more light will be seen with the _L polarization than with the || polarization (for which 
the image is unchanged). An observer of a highly magnetized NS will see a more complex image in which the same surface element can 
appear twice (or even more). 



multiple light ray solutions with the same v g (3 instead of 
!)• 

If the NS has a uniform surface temperature, its lumi- 
nosity will just be proportional to the surface area observed. 
This apparent surface area as a function of the magnetic in- 
clination i of the observer can be well fitted with the follow- 
ing functional form: 



4 

A' = a + J~] a„ cos"(2z). (8) 

One can also define an asymmetry parameter as the 
normalized difference between the area observed above the 
equator and that observed above the pole: 

= A'(i = 0°)-A'(i = 90°) (ai+a 3 ) , . 

A'(i = 0°) + A>(i = 90°) ' a {> 
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Figure 3. The fit parameters and the polar/equatorial asymme- 
try describing the image area for an isothermal NS. 



The value of the fit parameters and the total asymmetry are 
shown in figure ^. Evidently, one finds simple expressions for 
the fitting constants in the limit of weak fields: 



A « 0.095/q = 



S 



for Bo ^ 10 16 G and: (10) 



1.20 x 10 18 G 

ffl! « 0.89A, a 2 w -0.08.4, a 3 ^ 0.1U, a 4 « -0.044. 

With this fit, we can construct a light curve for a thermally 
uniform rotating NS. 



4.3 Average properties 

When a surface element is emitting black-body radiation 
into vacuum, we know from Stefan-Boltzmann's law that 
the emitted flux is simply F = uT 4 . Therefore, the intensity 
of each ray is I = B{T) = aT 4 /ir. This is however changed 
when the element is emitting the black-body photons into a 
medium with a different index of refraction. The reason is 
that in the presence of n / 1, there are on one hand more 
photon states for a given energy band and on the other, the 
group velocity which relates the density of states to the flux 
is smaller. For a nondispersive medium with an isotropic in- 
dex of refraction, the two factors are n 3 and 1/n respectively, 
implying that a unit surface will emit a flux of F = aT 4 n 2 , 
as if the surface area is effectively larger. We can therefore 
define an effective surface area as the equivalent area rele- 
vant for the Stefan-Boltzmann law of the emitted flux if the 
object would have emitted the photons into vacuum. That 
is: 



A n 



aT 4 



-dA. 



(11) 



where yl max is the effective surface area (and the notation 
will soon be clear). 



A careful analysis shows however that there there is a 
second relevant definition of an effective area. If an objects 
of surface area A and temperature T emits radiation into 
vacuum, then the total luminosity observed through a dis- 
tant surface surrounding the <in sterad around the object, 
will simply be: 



L = AoT 



(12) 



We can therefore define an effective surface area when the 
index of refraction is not unity as: 

L 

Tf 1 ' 



A 



(13) 



Under normal circumstances, we expect that the total flux 
emitted by a surface is the total flux observed at infinity. 
Namely, that A e s = j4 max . However, this need not be the 
case if some of the light trajectories leaving the surface in- 
tersect the surface again, or in other words, if there is photon 
trapping. Generally, we can expect that: 



A e ff < A ma x(= An 2 if n is isotropic). 



(14) 



There are two effective surface areas and both of them are 
interesting. The messy calculation of A max , the maximum 
theoretical effective surface area, is found in the appendix. 
Its results are compared with the observed effective surface 
area A e a calculated here. 

The effective area A c s of the star which is the total area 
as observered from infinity, is related to the apparent area A' 
seen at various inclinations through the following average: 



A cS = 4 



71-/2 



A'i 



i(i)di 



(15) 



The definitions are such that in the unlensed limit, the im- 
age has A' of tyR 2 while the surface areas A max and A e g are 
both A-nR 2 . The effective surface area can be directly related 
to the expansion of A' found in eq. |^ through its integration, 
giving that A oH /4 = ao+ai/3+7a 2 /15+9a 3 /35+107a4/315. 
It is however calculated directly in order to obtain high ac- 
curacy. 

A e ff is calculated using the ray tracing algorithm. For a 
given magnetic field -Bo and an observer inclination i above 
the magnetic equator, the apparent surface area A' is cal- 
culated. Since the surface brightness of the isothermal NS 
is the same in each direction in which light rays intersect 
the NS, one can achieve a high accuracy without the divi- 
sion of the screen into many pixels but instead use a more 
sophisticated algorithm. Instead, the accuracy is achieved 
through the accurate measurement of the apparent 'edge' 
of the NS at different angles relative to the center of the 
screen (through consecutive bisections) and then an integra- 
tion over angle gives an accurate area measurement - A'(i). 
This area can in turn be integrated to give A e s. 

The result for different magnetic field strengths can be 
seen in figure |i| which includes also the total emitted radia- 
tion from the surface A max as calculated in the appendix. 

Clearly from the figure, photons can only be trapped 
when the equatorial magnetic field is larger than 2.8 x 
10 16 G, corresponding to the point at which nontrivial dis- 
tortions of the apparent images of the NS take place. This 
also corresponds to the point where the contribution to n of 
higher order diagrams than the one loop approximation may 
be significant. Hence, this region should be taken cautiously. 
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Figure 4. The effective area of the NS when observed with the _l_ 
polarization. For equatorial fields larger than about 5x 10 16 G, the 
effective area has two definitions: the first is defined through the 
total emitted flux while the second is defined as the flux observed 
at infinity. The former can be larger than the latter when some 
photon trajectories are trapped. For small fields one finds a linear 
relation between the increase in the effective area and the field 
strength. 



For smaller fields, corresponding to the known field 
strengths of even the most highly magnetized PSRs, the ap- 
parent surface area in the _L polarization is: 

B 



An 



All 



5.0 x 10 15 G 



for B £ 10 iB G.(16) 



4.4 The light curve 



Generally, the axis of rotation and the magnetic field are not 
aligned. This implies that during the rotation of the star an 
observer will view it from a varying magnetic inclination. If 
the angular separation between the magnetic axis and the 
rotational axis is 7, then one can relate the magnetic incli- 
nation angle i to that of the inclination above the rotational 
equator i r and the rotational angle <f), which is the rotational 
phase between the last time the two axis coincided in the 
observer's meridional plane. The relation is: 



sin i = sin i r cos 7 + cos i r sin 7 cos < 



(17) 



(Greenstein & Hartke 1982). The best fitting expansion for 



the apparent area (and therefore the apparent brightness 
of an isothermal NS as well) was found to be a series of 
cos(2i). When expressed directly with 7, i r and <f>, one then 
finds that: 



A'(cf>) = ao + Va 



(l — 2 sin 2 i r cos 2 7— (18) 
2 cos 2 i r sin 2 7 cos 2 <f> + sin(2i r ) sin(27) cos 0)™ 



n=l 
2 • .2 



Since the odd terms in the expansion dominate the even 
terms (see eq. p"o| ), the expansion is predominantly an odd 
function of the parenthesized term and the temporal be- 
haviour (i.e., the maxima and minima) will be determined by 
this term as well. A simple analysis shows that for \i r \ > 7], 
there is one minimum and one maximum in one rotational 
cycle. However, for \i r \ < there are two equal peaks 
and two different minima in one rotational cycle, as can be 
seen in figure ^. If we neglect the contribution of the even 
terms, as we can for the weak field limit, the amplitude of 
the variations is just: 



A' max - A' min = aiFi(7,v) + 03*3(7, ir), 

where we have defined two auxiliary functions as: 



Fi(7,ir) = 



#5(7, *r) = 



2sin(2i r ) sin(27), 
2sin 2 (i r + 7), 



> |7| 
< 7 



(19) 



(20) 



i (3sin(2i r )sin(27) 

+ sin(6i r ) sin(67)) , 
(3 + 2cos(2(ir + 7)) 
+ cos(4(ir. + 7))) sin 2 (ir + 7), 



> M 



< 7 



The two functions incorporate the two regimes for which 
there are either one minimum or two minima in one rotation 
period. The difference between the depth of the two minima 
in the second case (when \i r \ < | — y | ) is the amplitude found if 
the \i r \ > I7I regime is used for the evaluation of F\ and F3. 
If ai is negative, the peaks become valleys and vice versa. 

A rotating isothermal neutron star observed in unpo- 
larized light will therefore have light variation of which the 
amplitude in the weak field limit is given by, 



AL 



4' 

- 1 in 



A' 

1 1 min 



A + A aV r 

(Fi(7,i r )+0.12F 3 (7,»r)) ( 



(21) 



B 



6.3 x 10 17 G 



The functions Fi and F3 are linear for small 7's and 
linear in the colatitude n/2 — - i r . For general values of 7 
and i r , they are of order unity. Namely, for small magnetic 
fields, the effect will fall as oc 7B 2 when viewed at random 
viewing angles and small separations between the axes. A 
few sample light curves are seen in figure ^. 

A more interesting effect however arises when studying 
the polarized light curves. In the unpolarized case, the effect 
stems from the asymmetry between the polar image and the 
equatorial image, an asymmetry that is never large, even for 
fields where n — 1 approaches unity. If we, however, look in 
polarized light, we have two images of the neutron star of 
which the area difference can be large for intermediate fields. 
Moreover, we shall soon see that irrespective of 7, the effect 
is always large and falls only as ~ B for small fields (as long 
as B <; -Bqed). Thus, it is expected to be more prominent. 

Using the polarization locking behaviour of the light 
rays (see §3.1), we know that photons leaving the surface 
with a || or _L polarization will keep their polarization until 
a large distance from the star at which point the local mag- 
netic field lines can be assumed to be totally aligned with the 
apparent direction of the magnetic axis (with the direction 
of the axis' projection onto the observer's plane of the sky). 
Since the magnetic axis is rotating, the two primary polar- 
ization directions are the projected direction of the axis onto 
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1.5 0.5 
t/Period 

Figure 5. A few sample light-curves of an isothermal NS having 
an equatorial field of 10 15 G, observed in unpolarized light, for 
different magnetic axis - rotational axis separations 7 and differ- 
ent rotational inclination angles i r . 



the plane of the sky and its perpendicular direction within 
this plane. 

If we work in a coordinate system aligned with the ro- 
tational 2-axis, and a j/-axis that is perpendicular to the 
plane containing the line of sight and the z-axis, then the 
observer's direction is: 



cos i r x + sin i r z 



(22) 



If we use the rotational phase <j> and the separation 7 be- 
tween the two axes, the direction of the magnetic axis is: 



rh = sin 7 cos <f> x + sin 7 sin <f> y + cos 7 z 



(23) 



Using this, we can calculate the cosine of the apparent angle 
zu that the magnetic axis makes with the y-axis (it is also 
the horizontal direction of the observer's screen). To do so, 
we project the magnetic axis to the plane of the sky and 
calculate its dot product with the observer's y-axis (which 
is also our coordinate system j/-axis): 



2 

cos zu 



((rh - (rh-6)6) ■ y) 2 
|rh — (rh • 6)6| 2 

(sin 7 sin < 



(24) 



:p(l,ir,4>) 



1 — (cos 7 sin i r + cos i r sin 7 cos (j> 

If we neglect the slight asymmetry between the apparent 
size when observed from the poles or from the equator, then 
the images observed at one or the other polarizations that 
are aligned with the magnetic axis are going to be fixed at 
A and A c s. If however our two observational polarization 
directions are fixed at for example the rotation axis and the 
perpendicular direction, then the polarized light-curve ob- 
served will vary because of the varying component observed 
from each polarization. 



The brightnesses of the two polarizations fixed to the 
magnetic axis are 

A 



L|| = 
L± = 



A + A 
A 



A + A 



-L w i (l - £p\ L for AA « A (25) 
iff 2 \ A J 

h L ~K 1 + ^r) L fOT * A « A ^ 



The two light curves at the two polarizations will therefore 
be: 



L x = 



Lj_ cos 2 w + iii sin 2 zu 



L 
2 



L AA 
2~ 



(2p(7,tr,^)-l). 



L2 = L± sin 2 zu + L« cos 2 zu 

« §-~^(2 P (7,ir,*)-l). 

where for small fields one has 
AA _ B 
~A~ ~ 5.0 x 10 16 G 



(27) 
(28) 

(29) 



This light curve normally has two peaks during one ro- 
tational cycle but the peaks do not have to be symmetric. 
The averages of each polarized light curve need not be equal, 
as can be seen in figure |(| However, irrespective of whether 
there is an asymmetry and irrespective of 7, one normally 
expects to see two different average amplitudes for the two 
different polarizations. 

Another interesting point is that even if the magnetic 
and rotational axes are aligned and there is no variation at 
all of the light curves during rotation, one still finds that the 
two polarizations will have a different amplitude: 



AL 



L2 — Li 



L± - Li, 



AA 



(30) 



5 A REALISTIC TEMPERATURE 
DISTRIBUTION 

Neutron stars are not in thermal equilibrium. They are born 
with T J> 10 9 K ((Shapiro fc Teukolsky 198^). After about 



one hundred years, their cores become isothermal and cool- 
ing proceeds quasistatically. Until 10 s — 10 6 yr have passed, 



the neutrinos dom i nate the heat flux from th e core ( Shapiro 
& Teukolsky 1984 |Heyl fe Hernquist 1998b| ). Later photon 



emission from the surface carries heat from the neutron star 
Regardless of the cooling process dominating at a particular 
time, the energy to supply cooling emission from the sur- 
face must traverse the thin insulating layer surrounding the 
neutron star core, the envelope. For warm and hot neutron 
stars the mo st important vehicle for the heat flow is e lectron 
conduction (Gudmundsson, Pethick & Epstein 1982). 

A strong magnetic field dramatically alters the conduc - 
tion properties of neutron star envelopes (Hernquist 1984). 
Electrons flow more easily along the field lines than across 
them, which causes th e flux distribu tion across the surface 
to be anisotropic (e.g. Schaaf 19908 



Schaaf 19901 



Hernquist 1998a, Heyl & Hernquist 1998b). For magnetic 



fields stronger than 10 14 G, analytic techniques can model 
the t hermal structure of the en velope with adequate preci- 
sion. Heyl & Hernquist (1998a) find that a simple prescrip- 



Heyl & 
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Figure 6. A few sample light-curves of an isothermal NS having 
an equatorial field of 10 15 G observed in polarized light, for differ- 
ent separations (7) between the magnetic and rotational axes and 
different rotational inclination angles i r . The two curves are the 
observed light curve in the two polarizations parallel and perpen- 
dicular to the rotational axis. (Solid line - when the ray's B field 
is parallel to the rotational axis, dashed - when it is perpendicular 
to the axis. 



tion describes the crust flux distribution for these highly 
magnetized neutron stars: 



F c oc B a4 (s- 



(31) 



where s is a unit vector normal to surface, and B and B are 
the magnitude and direction of the magnetic field. 

This indicates that in addition to the effect of the asym- 
metry between the polar image and the equatorial image 
that arises from the magnetic lensing, a large asymmetric 
contribution will arise from the flux gradient along the sur- 
face. Hence, it is very likely that the latter asymmetry will 
drown the unpolarized effect of § 4.4 . 



When calculating the emitted fluxes, one has to take 
into account that it is the flux driven through the crust 
which is fixed constant and not the surface temperature. In 
other words, a surface element will emit a total flux given 
by eq. |3l[ even if it is a more efficient black-body radiator. 

Due to the increased number of states in one polariza- 
tion (per unit energy) and the reduced group velocity, the 
black-body radiation law becomes: 



F x =aT 4 (l + ff (B,V))/2 and F n = aT 4 /2, 



(32) 



where the function g(B, ip) describes the increase in flux due 
to the magnetic field strength B that makes an angle tp with 
the surface normal. A fitting function for it in the limit of 
weak fields is found in the appendix. 

If the equilibrium temperature without the magnetic 
field would have been To, then it is evident that the new 



equilibrium temperature is going to be smaller in order to 
keep the same outgoing flux constant: 

T " = l+ 5 (zU)/2- (33) 

When we now ray-trace to find the resulting image, we 
have to take into account that the flux driven into the surface 
is a function of the magnetic field and that the temperature 
is a function of the field as well. Each ray coming from the 
observer's screen that intersects the surface will have, irre- 
spective of the polarization, a radiation intensity of: 



B 



cos 2 ip 



■Bpole/ 1 + <7(B,V0/2 7T 



-oTr 



pole 5 



(34) 



where T po i c is the polar effective black-body temperature. 
Excluding the index of refraction effect that enters through 
the denominator, the flux emerging from the poles will be 
larger than that of the equator. As a consequence, the image 
observed from the polar direction will be brighter than the 
image from equatorial inclinations. The effect of the index of 
refraction is to slightly cool the star. Since the image in the 
||-mode is unaffected by the nontrivial n, it will be fainter. 
On the other hand, the increased image surface area of the 
_L-mode will imply that it is brighter than the unlensed im- 
age, even though the lensed star is cooler. If averaged over 
the 4-7T sterad around the star, the sum of both images should 
yield the same total flux that it would have had without in- 
cluding of the refraction effects. 

To find the light curves, we find fitting formulae to the 
total image brightness as a function of magnetic inclination. 
They are accurate to roughly 1 part in 10 3 for very small 
fields. First, the we fit image brightness in the limit of B —> 
0. It is: 



L..(B-»-0) = Lx(B->0) = Lo &o 



4 



6„cos n (2i) (35) 



b = 0.5205, 61 = -0.1318, b 2 = 0.011, 
63 = -0.0026, and 64 = 0.0007. 

Lo o-Tp olc is the brightness that would have been observed 
if the NS would have had no magnetic field and would have 
had an isothermal temperature distribution with T po i e . 

Next, we fit for the corrections to the observed bright- 
nesses for each of the two polarizations: 



4 



L|| (B -> 0) + L 6||, + V&||,n cos n (2i) 



n=l 



B 
Bo" 



(36) 



6 ||0 = -0.7792 



5 II,3 



by,! = 0.2546, 6|| l2 
0.0001, and 611 4 = 0.0048, 



-0.0077, 




Ll = L ± (B 0) + L I 6 x ,o + ^6 

\ n=l 

&x,o = 0.4069, &x,i = 0.4960, b ± , 2 = 

b±, 3 = -0.0429, and b ±A = 0.1056. 

With these fits, we can construct any light-curve in the 
limit of small magnetic fields. The polarization brightnesses 
vary by more than 50% when observed at different latitudes, 
however, if we look at the ratio between the two, we single 
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Figure 7. The weak field ratio between the total flux observed 
in two polarizations as a function of the observer's magnetic in- 
clination i, normalized to /o- For small magnetic fields, the ratio 
Lj_/L^ is 1 plus the graphed function times /o. 

out the effect of the magnetic refraction. The ratio (for small 
magnetic fields) normalized to a magnetic field strength of 
Bq is depicted in figure 0. For 10 15 G, one finds that the 
ratio varies between roughly 1.015 to 1.03 when the mag- 
netic inclination is varied from the pole to the equator. The 
reason that the equatorial image brightness ratio is larger 
than the polar one is because the "over-the-horizon" effect 
of the magnetic refraction affects the equatorial image more 
because the luminous poles lie on the horizon - its inclusion 
is more important than the inclusion of the faint equatorial 
regions in the polar images. 

To construct the light-curves as a function of i r - the 
inclination above the rotational axis, 7 - the separation an- 
gle between the two axes, and <j> - the rotational phase of 
the star, one first needs to calculate for given angles the ex- 
pression for cos(2i), as is given in eq. hjJ. This allows the 
calculation of L« and L± through the use of eqs. [§6| and [Sri 
Then, using eq. M one can obtain the polarized light curves 
when observed in the two directions that are aligned with 
and perpendicular to the rotational axis of the star. A few 
sample polarization ratio light-curves are found in figure || 
One set of curves describes the total polarization, i.e. , when 
the measurement directions follow the || and _L polarizations, 
while the second set of curves is a more realistic case in which 
the observer's polarization planes are fixed and aligned with 
or perpendicular to the rotational axis. 



6 DISCUSSION 

The magnetic field surrounding a neutron star magnifies and 
distorts the image of its surface. If the surface field is every- 
where less than 3iv/a x Bqed ~ 5.7 x 10 16 G, the distortion 



Figure 8. Several sample polarizations light curves predicted for 
a 10 15 G neutron star. For different viewing inclinations i r above 
the rotation equator one will observe a changing ratio between 
the two polarizations. The ratio L±/L^ (solid line) is the flux ra- 
tio of the two polarizations while having the observer's polariza- 
tion planes follows the two polarizations states. The L1/L2 ratio 
(dashed line) is the ratio predicted if the observer's polarization 
directions are aligned with the rotational axis of the neutron star. 



is trivial; the neutron star appears slightly larger, and some 
light from the far side of the star contributes (the left side of 
figure ^). In this regime, only the image in the perpendicu- 
lar polarization is magnified. If each element of the neutron 
star's surface radiates as a blackbody, it will appear 3.2 % 
brighter in the perpendicular polarization than in the paral- 
lel polarization for the inferred surface field of SGR 1806-20. 
The flux difference between the two polarizations increases 
linearly with the strength of the field at the surface. 

Furthermore, as the neutron star rotates, the image in a 
fixed polarization consists of a varying admixture of the im- 
ages of the parallel and perpendicularly polarized photons. 
The phase of this oscillation at a given time depends weakly 
on the energy of the photons. The observed polarizations of 
higher energy photons reflect the polarization modes of the 
neutron star slightly later in its rotation period (eq. ^). 

In regions of the magnetosphere where the field strength 
exceeds 5.7 x 10 16 G, the index of refraction for perpendicu- 
larly polarized photons will exceed 3/2, and the relationship 
between the group and phase velocities of the light rays be- 
comes nontrivial. Regions of the star appear in the image 
several times with differing parities. The index of refraction 
is neither isotropic with respect to the direction of the pho- 
ton's propagation nor isotropic around the star; therefore, 
the optically dense magnetosphere produces a complicated 
image of the neutron star surface (the right side of figure |J). 

Although the full complexity of the images created by 
supercritical magnetic fields surrounding neutron stars will 
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be difficult to observe in the near future, several hallmarks 
of these effects can be observed with some of the next gen- 
eration of X-ray telescopes (e.g., Spectrum- X-Gamma), if 
sufficient precision can be achieved to observe the effects 
averaged over the two polarizations. If a particular type of 
emission is restricted to small regions of the magnetosphere 
of the neutron star {e.g. curvature and synchrotron emission 
may be restricted to the polar regions), this emission may 
be strongly magnified even in weak fields. Figure [] depicts a 
neutron star whose apparent area is about 50% larger than 
its actual area. The areas above 80° latitude are magnified 
by factor of up to 2.3 depending on inclination. Generically 
the polar regions will be magnified twice as much as the 
star on average since the polar field is twice as large. In the 
strong field limit, the complex topology of the images would 
be reflected in a complicated light curve for emission local- 
ized to the polar regions. Several pulses with varying parities 
could be visible from each polar region as the star rotates. 

Time-resolved spectropolarimetry of strongly magne- 
tized neutron stars would reap substantial rewards. All of 
the effects discussed in this paper affect only the perpen- 
dicular polarization mode. High energy photons traveling 
in the parallel mode will suffer photon splittin g and trans- 



fer their energy into th e perpendicular mode (Adler 1971 
Baring 1995, Heyl 199q ). This second process will affect the 



spectra of these sources especially for E ~ m e c . The lens- 
ing process discussed here affects even low energy photons 
(E <C m e c 2 ) - the effect will be even stronger at higher 
energies since the i ndex of refraction increases with photon 
energy (Adler 1971). Future generations of X-ray telescopes 



may be able to measure not only the energy but also the 
polarization of the photons that they detect (R. Rutledge, 
private communication) . 

Even in weak fields, the perpendicular polarization from 
the polar regions would subtend a larger solid angle in the 
image than the parallel polarization; consequently, its pulse 
would last longer. Moreover, the gross magnifications of the 
star and the augmented magnification of the polar regions 
would change the light curves in the perpendicular polariza- 
tion both in magnitude and character. For a uniform tem- 
perature blackbody, the perpendicular polarization would 
be brighter than the parallel polarization and exhibit slight 
variability. The parallel polarization would be what one 
would expect for a blackbody at the appropriate temper- 
ature, modulo the effects of photon splitting. If the star has 
a nonuniform temperature distribution with the poles hot- 
ter than the equatorial regions, the lensing will increase the 
variability observed in the perpendicular polarization, and 
also increase the total flux detected in that mode. 

Several additional effects could quantitatively change 
the results, two of which are gravitational lensing and the 
presence of higher multipole fields. Although gravitational 
lensing does not compete directly with magnetic lensing in 
the sense that its effect is independent of polarization (and 
thus creates no polarization signal) , it can slightly smear the 
signal by adding over-the-horizon contributions. 

The effects of higher multipoles can on the other hand 
be potentially more important. Multipoles that are some- 
what higher than the dipole will be directly observable as 
they will contribute higher Fourier components to the po- 
larization light curve (e.g. fig. 8). Under favourable con- 
ditions, it could allow the measurement of each multi- 



pole separately. Very high multipoles will not have a vary- 
ing signal since the observation of a whole hemisphere 
will average their variations. Nevertheless, since an average 
field strength does exist, a constant polarization signal of 
L±/L\\ ~ 1 + min{B mu iti P oie/10 ir G, AR/R} will be intro- 
duced. AR is the typical height above the surface on which 
the high multipole field decays. It bounds the maximal L± 
since the increase of the apparent area cannot be larger than 
the cross-section of the region where the index of refraction 
is not unity. Although multipoles can complicate the anal- 
ysis, they could in principle be detected if they were more 
important than the dipole field. 



7 CONCLUSIONS 

We have studied how photons travel through the magneto- 
spheres of strongly magnetized neutron stars. For surface 
fields less than 3n/a x Bqed ~ 5.7 x 10 16 G, the mag- 
netosphere distorts the image of the neutron star surface 
minimally. The total area of the image in the perpendicu- 
lar polarization is mag nified by 5/(5.0 x 10 14 G) percent, 
and some light from the far side of the star contributes. For 
fields stronger than 5 x 10 G, the distortion of the image 
becomes more complicated. Regions of the stellar surface 
appear several times with varying parity in the i mage 



Unlike gravitational lensing of neutron stars ( Page 1995 , 
Heyl & Hernquist 1998a), magnetic lensing only affects the 



perpendicular mode and tends to increase the variability of 
the source; consequently, by measuring how the polarized 
flux varies as the magnetar rotates, one could observe the 
lensing described in this paper and estimate the strength 
and inclination of the surface magnetic field in a manner 
independent of the details of the dipolar spin down. 



APPENDIX A: THE MAXIMAL EFFECTIVE 
SURFACE AREA OF A MAGNETAR 



In §4.3 we have seen that there are two definitions for the 
effective area of a magnetar - the area required to emit the 
same black body luminosity by an object without a mag- 
netic field. The first is the effective area to an observer at 
infinity (using the total luminosity observed there) while the 
other is the effective area to the object itself (using the to- 
tal emitted luminosity). The definitions are not the same 
because some of the emitted photons can be reabsorbed by 
the star thro ugh trapped photon trajectories. 
In 



4.3 



we have calculated the effective area seen at 
infinity. To see whether photon trapping exists, one has to 
calculate the total luminosity emitted by the star. 

The effective area increase of a surface element is given 
by the relative increase in the black-body flux. Since the 
index of refraction n is not a function of energy, the flux 
increase is simply an angular integration: 



da_ 

da 



L 

f 



df2j.n(k) 



s) 



(38) 



v g (k)-s>0 



The integral is over the half sphere for which the group 
velocity v g points upwards, away from the surface (of which 
the normal is s). The n(k) 3 factor is the increase of the 
number of possible photon states for a given energy interval 
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Figure 9. The fc-space over which the flux integration should 
be carried is the region in which v s • s > — i.e. , above the 
hatched line. Each point in the integration region (e.g. 1,2) has 
however a counterpart outside the region (l',2') for which the 
group velocity is opposite. If we integrate |v s -s| over all space we 
get twice the required answer without the complication arising 
from the boundary. 



at a direction flk. The factor v g (k) • s is the contribution to 
the flux perpendicular to the surface from these states. For 
an isotropic nondispersive medium, one immediately finds, 
as one should, that /'// — ■ 

The integral boundaries might appear simple but they 
are in fact very complicated because v g is generally not in 
the direction of k. We can however simplify the boundaries 
through the use of a symmetry property of n(k). Since n 
is parity conserving: n(— k) = n(k) and since v g is a vec- 
tor that is related to n through a gradient, it will satisfy: 
v g (— k) = — v g (k). If we now observe figure H describing the 
vertical slice of the integration regions through the plane 
containing B and s, then it is clear that each point within 
the proper integration region (e.g. point 1 or 2) has an equiv- 
alent counterpart with an opposite group velocity outside 
the integration region (e.g. point 1' or 2'). Thus, to get the 
required integral, we can either integrate using the cum- 
bersome boundaries or integrate the absolute value of the 
integrand of the whole region and divide by 2. We choose to 
do the latter in the polar coordinates (r, 9, <p) aligned with B 
such that the y axis is perpendicular to the plane containing 
both B and s. 

In our coordinate system, we have that the radial unit 
vector is: 



f = sin 9 cos <p x + sin 9 sin <p y + cos 9 z. 



(39) 



The unit vector in the direction of 9 is needed for the group 
velocity, it is: 



= cos 9 cos 4> x + cos 9 sin (py ~ sin 9z. 



(40) 



Last, we need the vector in the direction normal to the sur- 
face element. If the angle between the magnetic field direc- 
tion and the normal to the surface is ip, then this vector 
is: 

s = — sin ip x + cos ip z- (41) 
The group velocity is related to uj through: 



^_ , . . ~ , 1 / duo 

V ^^={dk) e r -k{d9j k 



(42) 



If the index of refraction is in the form n - 
then we find: 

1 2 f(B) sin 9 cos 9 



1 + f(B) sin 2 6, 



3 l 

n v.. 



1 + f(B) sin 2 9 (1 + f(B) sin 2 9) 2 
The integrand is therefore: 

(l + /(B) sin 2 9) 2 x 

(cos ip cos 9 — sin ip sin 9 cos <p) 
+2f(B) sin 9 cos 6 (l + f(B) sin 2 0) 
(sin ip cos 9 cos <p + sm 9 cos ip) . 



(43) 



(44) 



Next, we need to integrate over the NS surface with 
its varying field strength and varying directions. Since we 
assume it to be a dipole, we have that if the equatorial mag- 
nitude of the magnetic field is B eq then on the surface we 
have: 



B = B cq (3(s ■ rh)s — rh) . 



(45) 



We thus find that the magnitude of B and its angle from 
the local normal are: 



B = 
cos ip — 



B s 



m) 2 + lB cq 
2 cos a 
\/3 cos 2 a + 1 



\J 3 cos 2 a + lB a 



(46) 
(47) 



We further use that f(B) is linear in B in the interesting 
regime, implying that f(B) = f(B cq )B = /o\/3 cos 2 a + 1. 

The growth of the effective area of the whole star can 
now be written as: 



da 4-7T / da 



(48) 



where J dQ, a implies an integration over the surface colati- 
tude and longitude (a,ip). The integration over <p is trivial 
and we have: 

1 ^ — 1 I.2TT 

d(cosa) / d(cos9) / d<p 



1 

4tt 



(l + f(B) sin 2 (cos ip cos 9 — sin ip sin 9 cos <p) 

+2 f{B) sin 9 cos 9 (l + f(B) sin 2 0) x 

(sin ip cos 9 cos <p + sin 9 cos ip) | , (49) 

where we explicitly have to express f(B) as /o V 3 cos 2 a + 1 
and cos ip as 2 cos a / V3cos 2 a + 1. The integration is non- 
trivial (mainly because of the absolute sign) but is however 
straight forward if done numerically. The method used was 
a Monte Carlo integration which easily achieved accuracies 
better than 10 -3 for several dozen magnetic field strengths 
in roughly an hour of an HP workstation CPU time. 

Another useful computation is the direct evaluation and 
the construction of a fitting formula for the local increase 
in the black-body flux as given by eq. The integral is 
generally a nonlinear function of both the magnetic field 
strength B and the angle ip between the magnetic field, B, 
and the surface normal, s. Thus, a simple fit can only be 
found in the limit of weak fields (B < B w 1.13 x 10 17 G), 
it is: 

, „ , % da' 
I + 9(B,1P) = - 



1 + (l + 0.067 cost/) + 1.6024 cos 2 ip 



-0.672 cos 3 ip) (B/B ) 



(50) 
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Its relative accuracy is better than 10 3 , which is the inte- 
gration accuracy. 
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